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By using of a special reduction way of density matrices, in this Letter we find the entanglement 
between two bunches of particles, its measure can be represented by the entanglement of formation. 
PACC numbers: 03.65Ud, 03.67-a, 03.65Bz, 03.67Hk. 
In this Letter, we discuss a special way of reduction of density matrices, and prove that in an arbitrary multipartite 
qubit state there is the new kind of entanglement, i.e. the entanglement between two bunches of particles, independent 
what happens to the remaining particles, and its measure can be represented by the entanglement of formation[l] Ef. 
, Some examples are discussed. 

It is known that the problems of the measure of the entanglement of bipartite qubit (pure- and mixed-)states are 
solved, e.g. see [1,2]. However, the description of entanglement of multipartite qubit states is a formidable task as 
£NJ ■ yet. Recently many new results have been obtained, one way in which is to describe the properties of entanglement of 
,_h j multipartite qubit states by using of the bipartite reduced density matrices and of the entanglement of formation Ef. 
For instance, in [3] the entanglement between two particles, independent what happens to the remaining particles, is 
described by a bipartite reduced operator, and its measure is represented by Ef. Some more new results concerned, 
GN ! see [4,5]. 

For the spin particle M^(k = 1, 2, • • • , N, N > 3), we simply write f Mk >=| >fc and J, Mk >=| 1 >k, or in union 
by | i >k {i — and 1). | i >fcspans the Hilbert space H^. The main points of the above ways are as follows: If p is a 
density matrix acting upon Hi <g) H2 ® • • • <S> -ffjVj let the pki be the bipartite reduced density matrix defined by 

*Q i pu=tr x x <j>), (l<k<l<N) (1) 

l... It— l —N 

j — ' where the tr denotes the trace of a matrix, and symbol x denotes the deletion operator. Since pki is a bipartite 
density matrix, we can use the entanglement of formation Ef. Therefore once use pki and Ef [pki] to describe the 
C*") ■ entanglement status between two particles | i and | i independent what happens to the remaining particles. 

About the above kind of ways by using of reduced density matrices pu and Ef [pm], we need to consider at least 
the following two problems: First, why only consider we two particles, but not more particles? In fact, for the case 
O ,' of entanglement among more particles the problem backs again to the original status, i.e. we need to handle other 
multipartite qubit state, the above way will runs up against difficulties. For instance, when N > 5, and we need to 
consider how to describe the entanglement among four particles | i >i, | i >a, | i >3 and | i >4, independent what 
happens to the remaining particles | i >5, • • • , | i >n, then we use the reduced matrix P1234 = f?*5...jv(p) 5 however P1234 
is not a bipartite qubit state, we cannot use Ef. Although we still can write the set {Ef [pi2\ , Ef [P13] , • • • , Ef [P34]}, 
it cannot show more contents of above 'entanglement among the particles | i >i, | i >2, | i >3 and | i >i, but the 
latter must contain other more contents. How are we to surmount this difficulty? Secondly, for some important 
multipartite qubit entangled states, say |GHZ W >= 4- (| 00 • • ■ > + | 11 • ■ • 1 >) , pghz n =|GHZ w xGHZjv |, all 
Pki are disentangled, i.e. all Ef [pu] — 0, this is somewhat making one puzzled: Is there some possible entanglement 
shared between two parts of system |GHZat > with non-zero Ef? In this Letter, we suggest a way that ones should 
consider some special ways of reduction of density matrices, then we prove that there is a new kind of entanglement 
in multipartite qubit systems, i.e. the entanglement between two bunches of particles, its measure can be represented 
by Ef. These results bring to light the more properties of multipartite qubit entangled states. 

In the first place, we discuss the simplest case, i.e. the tripartite qubit entangled states. The general form of a ^ in 
H a (g)H b (g) H c is as * = J2i,j,k,=o,i I i >a ® | j >b ® | k > c ( cy* G C). We define H a/bc , -ff a/fe v , respectively by 

H a/bc = < * |: The form of * is as * a/bc = ^ c lkk \ i >a <S> \ k > b | k > c > (2) 

[ i,fc=0,l J 

H ,. v = I * |: The form of * is as * ,. v = V c ife( i_ fe) | % > a ® | k > b ® | 1 - k > c 

a/bc a/bc i. — ^ v ' 1 

i,fc=0,l 

H a / bc , H /t v are two 4-dimensional subspaces orthogonal to each other, we have the direct sum decomposition 
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a/bc 



H a ®H b ®H c = H a/bc 8 H a/b . (3) 
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Now wo make two formal bases I i > x and I i > y , since H a / bc and H v both are 4-dimensional spaces, we take the 

' a/bc 

1-1 correspondences as 

I i > a <8> | k > b ®k > c ^| i > a ® I ^ >x and | i > a ® | fc >& ® | 1 — k > c ^| i > a ® | fc > y (4) 
then we have the following isomorphisms 

Ha/bc ~ H a ® H x , H w ^ H a ® H y , H a ®H b ®H c ^H a ®H x ®H a ® H y (5) 

From Eq.(5), any quantum state \E< can be expressed in only one form as a sum of two orthogonal states ^ a /bc 
and * , b v, especially, we can explain ^ a /bc and * , b v as follows: ^ a /bc is the 'wave function describing two bunches 

(| i > a ) and (\ k > b (g> \ k > c )', where in bunch (| k >b ® | k > c ) the spin-directions of particles b and c always are 
the same, hence their behavior of spin can be regarded, as a whole, like to a single spin particle. W is the 'wave 

function describing two bunches (| i > a ) and (| k >b ® | 1 — k > c )\ where in bunch (| k >b ® | 1 — k > c ) the spin 
directions of particles b and c always are contrary, hence their behavior also can be regarded, as a whole, like to other 
single spin particle. 

The projection from H a ® H b ® H c to H a (g> H x and H a ® H y , respectively, are two left multiplication operators as 
Pa/bc = I i >a I k >x<c k\® < b k\® < a i\:^ — > V ax = P a/bc (*) = Y c ikk \ i>a®\k> x (6) 

i,fe=0,l i,fe=0,l 

p a/b -= \i> a ®\k> x < c l-k\® < 6 fc | ®<„ i |: * — > = P Q/fe v (*) = Y c tk(i-k) \ i>a®\k> y 

i,fe=0,l i,fe=0,l 

And the interior mappings J ax : H a <g> P^ — ► P a (g) H b ® H c and / aj) : P a <g> Pj, — ► P a (& H b (g> H c are two left 
multiplication operators as 

/ aa; = ^ |t> ®|A;>6®|A;> c < x A;|®<„*|:*ax — ► ^ox (*«r) = *o/6c- (7) 

j,fe=0,l 

J 01/ ee ^ I % > a <g> I A; >6 <g> | 1 - k > c < x fc|®<„i|:* ow — ► I ax (tf oa: ) = * a/fc v. 

j,fe=0,l 

Obviously, P a / bc ° Prr and o J ai/ , respectively, are the identical mappings upon H x and H y . 

Suppose that T is an linear operator, 

T: H a ® H b ® H c — ► H a ® H b ® H c , $ = T (*) (8) 

then we can obtain the induced mappings T ax and from T, 

T ax : H a ®H x — > H a ® H x , ^ a x — ► ^ax = T ax (* ox ) = P Q / 6c oToI ax 

T ay : P a ® — ► H a ®H y , t>ay —> ^ax = T ay ay ) = P a/fc v ° T o I ay (* oy ) (9) 

We take especially an interest in the case of that T is a ( pure or mixed) density operator p on H a ® H b ® H c , in this 
case the results obtained are two bipartite operators T ax d and T ay d as 

• ^arcd * ^arc = (^ax) 3 ^aa: — Pa/bc ° P ° lax 

T ay : ^ayd — = T „ (*„„) , T ay EE P^v opol^ (10) 

By using of Eqs.(6), (7) and (9), we find the entries of T ax and T ay , respectively, are 

\^ax\ij t kl = [P\ijj,kll ' \^ a 'y\i],kl = \P\ij(l-j),kl(l-l) ' (*>.?> M = 0,1) (11) 
where [p] i j m Hs are the entries of density matrix p. By normalization, we can write 

"f ax = VaxPax, Vax = Mrss.rss ' = ~ T"a2:- = rj ay Pay, Vay — [p]rs(l-s).rs(l-s) ' Pay = ~ ay 

r, S =0,l ^ aa; r,s=0,l % « 

(12) 
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Since p is a density matrix, from Eqs.(10),(ll) and (12) we can directly verify that p ax and p ay both are bipartite 
density operators. Obviously, p ax describes the status of entanglement between bunches (| i > a ) and (| k >b <8> | k > c ) , 
and p a y describes the status of entanglement between bunches bunches (| i > a ) and (| k >b <8> | 1 — k > c ) . In addition, 
there is the relation rj ax + rj ay = 1. This means that we can consider the operator p( a . oc ) defined by 



P(a,bc) 



^ ax "i" f ay IJaxPax ~t~ VayPay 

Wooo.ooo + Wooi, ooi Wooo,on + Wooi, oio Wooo,ioo + Wool, 101 Wooo,m + Wooi, 110 

Won, ooo + Woio,ooi Won, on + Woio.oio Won,ioo + Woio,ioi Won, m + Woio,no 

Wl00,000 + Wl01,001 Wl00,011 + Wl01,010 Wl00,100 + Wl01,101 Wl00,lll + Wl01,110 

. Wm, ooo + Wno, ooi Wm, on + Wno,ooi Wm,ioo + Wno,ioi Win,ni + Wno,no 



(13) 



then p( a h c ) can be taken as a bipartite qubit mixed-state, which describes the status of entanglement 
between two bunches of particles (a) and (b, c). In addition, the procedure in accordance with the rules 
in Eqs.(ll), (12) and (13), in fact, is a special reduction of density matrices. 

Similarly, we take 



H 



b/ca 



^ I. The form of ^ is as ^ 



bj ca 



E 

i,fe=0,l 



Ckik I k > a ® I i >b ® I k > c 



H 



b/ca 



I. The form of * is as * 



b/ca 



i,fc=0,l 



i >b ® I k > c 



(14) 



Correspoces | k > a ® | i > b ® | k > c ^| i > b <Z> \ k > x =^ 
I l-k> a ®\i>b®\k > c ^\ i>b®\k> l 



H 



bj ca 

H, 



Hh <8> Hrr 



b/ca 



H h <g> H, 



b/ca 



and similarly construct the projections P b / ca ,P b ^v, the interior mappings Ib/ ca , ^ h/f v and the induced mappings 
Tfcx, Tf, y , • • • ,etc They lead to 

P(fc,ca) = TTfcr + T by , [p(b,ca) 



w, 



Wri- 



P(b,ca) 



ij,kl - VSjij,M "T m(l-j)ij,(l-l)kl 

Pooo.ooo + Pioo,ioo Poooaoi + Pioo,ooi Pooo,oio + Pioo,no Poooan + Pioo,on 

PlOl, 000 + P001, 100 PlOl, 101 + P001, 001 PlOl, 010 + Pooi, 110 PlOl, 111 + Pooi, 011 

Poio.ooo + Pno,ioo Poioaoi + Pno,ooi Poio,oio + Pno,no Poioan + Pno,on 

Pill, 000 + P011, 100 Pill, 101 + P011, 001 Pill, 010 + Pon, 110 Pill, 111 + Pon, 011 - 



(15) 



P(b.ca) is a bipartite density matrix which describes the status of entanglement shared between two bunches of particles 
(6) and (c, a). Similarly, we can yet obtain 



P(c,ab) 



T 

J- c 



- cy, 



[P(c.ab)] 



ij,kl 



[Pimm + Wj(i-j)i,i(i-ofc 



(16) 



Notice that although we can yet write P( a b,c), 1 ■ ■ i there are repeats, e.g. P( a b,c) = P(c,ab), 1 ■ ■ i etc., there only are three 
independent P(.,..), i.e. p( a ^ c ), P(6, ca ) and p( c . ab )- 

Since P( a ,bc)i P(fc.ca) and P( c ,ab) all are bipartite density matrix, we naturally use Ef to represent their entanglement 
measure. For a given p this measure Ef can be concretely calculated by using of the so-called 'concurrence '[6, 7]. For 
instance, for P( a ,bc) defined as in Eq.(13) 



Ef[p(a,bc)]= h (l + \Vl-C^ 



(17) 



where h is the binary entropy function h (x) = — xlog 2 x — (1 — x) log 2 (1 — x) , the concurrence C is determined by 

C = max{0,-Ai,-A2,-A 3 ,-A 4 ,} (18) 



where \ L are the eigenvalues, in decreasing order, of the Hermitian matrix R = \J ^P( a ,bc~ P(a,bc) \/P(aM)> P(a.bc) z 

-i~ 

1 

P(b.ca) and P( c ^ab)- Therefore, we can obtain the complete set of measures of entanglement shared between every pair 



(ct2 ® oi) (p(a,bc))* (°2 ® CT2) , (72 is the Pauli matrix 



, the star is the complex conjugation.. Similarly, for 
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of two bunches, i.e. {Ef [p( a ,bc)] ,Ef [p(b.ca)] ,Ef [p( c , &)]}> it describes some character of the entanglement status 
among three particles a, 6, and c. 

Now, we return to the problems mentioned in the start of this Letter. In the first place, if p is a density operator 
upon H\ <g> H 2 <8> • • • <8> Hn (N > 3), then tr x x x (p) (1 < j < fc < Z < N) is a tripartite density operator 

l—j—k—l—N 

acting upon Hi ® Hk <8> therefore the status of the entanglement between two bunches of particles 
(Mj) and (Mfc,M/), independent what happens to the remaining particles, can be described by state 



tr 



■N 



(p) 



. In the following we simply write 



UM) 



P(j,ki) 



tr 



■j — k— I ---N 



(p) 



(19) 



-I UM) 



Here we must stress that pu t ki) ls a special reduced matrix of p by two reduction procedures in succession: 
The first is the ordinary reduction, the second is in accordance with the special rules as in Eqs.(ll), 
(12) and (13). The measure Ef [p(j,w)] can be calculated as in Eqs.(17) and (18). Similarly, for P(k.ij) and P(i.jk)- For 
three particles | i >j, \ i >k, \ i >z,the set {Ef [p(j,w)] , Ef [p(k.ij)\ , Ef [p(i.jk)] } completely describes all entanglement 
between every pair consisting of a bunch containing single particle and a bunch containing two particles . Obviously, 
this shows a character of the entanglement among three particles | i >j, \ i >k, \ i >i in the multipartite qubit state 
p, independent what happens to the remaining particles. This cannot be obtained only by using of the ordinary 
reduction as in Eq.(l) and Ef. 

The generalization of more high dimensional bunches is straightforward, e.g. we can obtain P(i2,34)(^ > 4) from 
Pi2/34>P v ,p v,p v v, and obtain pn 234-1 from pi/234,p v ,p v,,p w (notice that, in fact, p v v = 

^' r 12/34 12/34 12/34 r ' ' r l/234 r l/234 1/234 1/234 

/9 i/234' e ^ c ")' Similarly, P(i3,24) 5 P(i4,23); ' ' "j etc At last, when N particles are divided into two bunches {ri, • • • ,r m } 
and {si, ■ ■ ■ , s n } , where 1 < n < r 2 < • • • < r m < N, 1 < si < s 2 < • • • < s n < N, {n, • • • , r m } n {si, • • • , s„} = 
and {ri, • • • , r m } U {si, • • • , s n } = {1, 2, • • • , TV} , then we obtain P({ ri ,---,r m },{si, •••,s„})- I n addition, the in Eq.(l) 
obviously is a special case of two 'bunches' containing only a single particle, or use our symbol, pki = P(k,i)- 

By the above ways, let (ki) m = {k\, • • • , k m } and = {h, • • • , l n } both be two subsets of {1, 2, • • • , N} , where 
m + n < N, 1 < ki < k 2 < ■ ■ ■ < k m < N, 1 < h < l 2 < ■ ■ ■ < l n < N, and {ki, • • • , k m } n {h, ■ ■ ■ , /„} = 0, then 
P((ki) (lj) ) is an entangled state between two bunches ({| ik t >k t , • • • , \ ik m >k m } , {| jh >i 1 ,--- , \ jl n >«„})• The set 

of all possible Ef p^ k .^ (;.) ^ (notice that there are repeats in jp((fc.) (;.) )|)j is yet a description of character 

of the entanglement among the m + n particles | ik 1 >fc x , ■ ■ ■ , | ik m >k m , | jh >h , • • • , | ji n >;„ in the N-partite qubit 
state p, independent what happens to the remaining N — m — n particles. 

Secondly, as a special example we consider the GHZ state pghz n ^|GHZat xGHZjv | (N > 3). By using of the 
above (h) m and {lj) ,we have the following results 



is discntanled, E 



p ((k,) m , (h) n ) is maximall y entanled, Ef [p( (fci)mi (; . )n ) 



0, for 2 < m + n < N 

1, for m + n = N 



(20) 



The proof only is a straightforward calculation by Eqs.(17) and (18). This result shows fully the character of pghz n1 
i.e. only when all N particles are divided into two parts (every particle must be in one and only one of them), the 
entanglement between this two parts does not vanish, and it is maximal. Therefore in view of this, for N> 3 the result 
that all Ef [{pGHZ N )ki\ = is not at all surprising. 

Other interesting example is {w is a given integer, 1 < w < N) 



V2 



0>i 



> t 



1 > 



w+1 



1 >N + I 1 >1 



1 >w ® | > w + i 



it like to the Bell state (p+ = ± (| > a O | 1 > b + \ 1 > a ® \ > b ) . For B 
verified that 



+ =1 

(N,w) "I 



ftw)«*,.. M «u---M) " entangled ' Ef 



M o >jv) 

(21) 
it is easily 



.!»}). 
M). 



= 0, for 2 < to + n < N 

> 0, for to + n = N (22) 
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More generally, if N < M and (rk) N = {ri, • • • , rjy} is a subset containing N elements in the set {1, 2, • • • , M} , we 
can construct </>t. , > and 5t. , x , according to Eq.(21) for the set {I it > rt } (A; = 1, N) of N particles. We 

™) (.( r fc)N > W ) Til 



define the state 



and the mixed state 



$ {M ,N, W ) >=\ <f>( (rh)ifi w ) > ® I ' ' ' >rest (23) 



B (M,N,w) = X (rk) N I *(M,JV,to) >< ®{M,N,w) I ( 24 ) 

all possible (r k ) N C{l,2,---,M} 

where the real numbers X( rk ) N obey 0<X( rk ) N < 1 and x (rk) N ~ 1> then Eq.(22) still holds 

all possible (r k ) N C{1,2,---,M} 

for Bf M Nw y The action of Bf M N , is somewhat like to an 'entanglement molecule' [3] . 

Sum up, in a multipartite qubit state p there is a new kind of entanglement, i.e. the entanglement between 
two bunches of particles, independent what happens to the remaining particles, which can be described by the 
special bipartite reduced density operators P(( fe .) ^ ), and the measure of entanglement can be represented by 

E f [P((M) m , (i,) B ) 
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